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INVESTIGATION OF SCENE IDENTIFICATION ALGORITHMS 


FOR RADIATION BUDGET MEASUREMENTS 
by 

F. J. Diekmann 

George Washington University, Hampton, Virginia 

ARSTRACT 

The computation of Earth radiation budget from satellite measurements 
requires the identification of the scene in order to select spectral factors 
and bidirectional models. In order to evaluate the effects of scene identi- 
fication errors on the resulting radiant exitances, scene identification 
probability matrices, which quote the probability that a scene which is in 
truth type j is identified as type i, are needed. This paper presents 
empirical evaluations of these matrices by intercomparing preliminary results 
from the ERBE scene Identification algorithm with results using AVHRR data. 

A scene identification procedure is developed for AVHRR SW and LW data 
by using two radiative transfer models. These AVHRR GAC pixels are then 
attached to corresponding ERBE pixels and the results are sorted into scene 
identification probability matrices. These scene intercomparisons show 
that there generally is a higher tendency for underestimation of cloudiness 
over ocean at high cloud amounts, e.g. mostly cloudy instead of overcast 
partly cloudy instead of mostly cloudy, for the ERBE relative to the AVHRR 
results. Reasons for this are explained. 

Preliminary estimates of the errors of exitances due to scene misidenti- 
fi cation demonstrates the high dependency on the probability matrices. 

While the longwave error can generally be neglected the shortwave deviations 
have reached maximum values of more than 12 % of the respective exitances. 
Numerical results in this paper must be upgraded with the final ERBE data 
when they become available. 
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I. INTRODUCTION 

The Earth Radiation Budget Experiment (ERBE) is a major effort in the 
understanding of the Earth's climate and its variations (Barkstrom and 
Smith, 1986). With ERBE, as in any satellite experiment for study of Earth 
energy budget at medium spatial scales, a scanning radiometer measures the 
radiance in one direction due to the upwelling refle'cted solar (shortwave) 
radiation or the Earth emitted (longwave) radiation. In order to compute 
the radiant exitance in either of these ranges, it is necessary to account 
for the isotropy of the radiation, which is described hy a limb-darkening 
function for Earth emitted radiation and a bidirectional function for 
reflected solar radiation. This approach was used for the Nimbus ?. and 
Nimbus 3 Medium Resolution Infrared Radiometer (Raschke and Bandeen, 1970; 
Raschke et al., 1973), the Nimbus 6 and 7 Earth Radiation Budget (W. L. 
Smith et al., 1978; Jacobowitz et al., 1984), and currently the ERBE (G. L. 
Smith et al., 1986). These functions depend upon the directions of the 
incoming and exiting rays, and also upon the scene, by which we mean the 
Earth's surface and the intervening atmosphere. The scenes are classified 
into a small number of scene types. The computation of radiant exitance 
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from the radiance measurement thus requires first the identification of 
the scene as one of these types, and the selection of the corresponding 
bidirectional and limb-darkening function. Due to the importance of these 
functions, their experimental determination as was done by Taylor and 
Stowe (1984) was one of the major objectives of the ERB experiment on the 
Nimbus 6 and 7 spacecraft. Furthermore, much computation has been done to 
define these functions theoretically, e.g. as by Stuhlmann et al . (1985). 

There is a second reason for the need for scene identification in the 
determination of Earth radiation budget. A scanning radiometer for Earth 
radiation budget ideally measures the broadband radiance in the shortwave 
and in the longwave ranges. However, the optics of any real radiometer 
will not have the ideal spectral response across the spectrum to as sharp 
cutoff, but will have a more complex spectral response. For example, the 
spectral responses of the three channels of the ERBE scanning radiometer 
are rather complicated functions of wavelength, as shown in Figure 1. The 
spectral responses of the shortwave and longwave channels of the ERB scanning 
radiometers aboard the Nimbus 6 and 7 spacecraft are essentially the same 
as these of the ERBE. For the ERBE data analysis, the spectral responses 
of the channels are accounted for by identifying the scene, which in turn 
provides an approximation of the incoming spectrum. From this information, 
the broadband radiances can be accurately computed (Avis et al., 1984). 

Because of the importance of scene identification in the data analysis 
process, it is necessary to understand the algorithm by which the scene is 
identified. This paper presents a study of the ERBE scene identification 
algorithm and helps to define its strengths and weaknesses. Similar con- 
siderations may be expected for any radiation budget measurements by sat- 
ellite. The approach in this study is to use radiative transfer theory to 
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establish a quantitative understanding of the various parameters of the 
problem. These theoretical results are compared with satellite data as a 
test of this understanding with real data. 

The performance of a scene identification algorithm is measured object- 
ively by the probability that a scene which is type i is computed to be 
type j. These probabilities for all combinations of i and j form a matrix. 
This scene identification matrix is shown by Smith et al . ( 1984 ) to be 
important in the calculation of error propagation and generation in the 
analysis of Earth radiation budget data. The present paper provides scene 
identification probability matrices on the basis of the comparison of ERRE 
derived scenes compared with the scenes as derived from AVHRR data. Because 
the scene identifications from the AVHRR are themselves subject to error, 
these probabilities are not with respect to absolute truth, but with respect 
to agreement of the two algorithms. Thus, as with any experimentally based 
quantities, these matrices are approximations. However, they are the first 
available. As such, they permit the computation for the first time of the 
errors in shortwave and longwave radiant exitances which are due to scene 
identification errors as measured by these scene identification probability 
matrices. Preliminary results for these radiant exitance errors are pre- 
sented. 

Because they are based on preliminary ERRE data tapes which are in the 
process of being revised, all results in this paper are preliminary. Indeed, 
a major reason for the work is to provide information concerning the per- 
formance of the ERBE scene identification algorithm so as to permit improve- 
ments to the algorithm prior to finalizing it for production of valid ERBE 
data. For this reason, the numerical results in this paper must be upgraded 
with the final ERBE data. 
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II. APPROACH 

We begin with a review of the ERBE scene identification algorithm and 
its role in data analysis. The basic requirements for a simulation study 
of the algorithm are thus defined, and the approach is designed. In order 
to study the application of the algorithm with real data, AVHRR data are 
used as "truth" against which the performance of the ERRE is compared. 

A. Description of ERRE scene identification algorithm 

The data flow for analysis of scanning radiometer data is set up as 
shown in Figure 2. Using "global" spectral factors, approximate broadband 
shortwave and longwave radiances L at the spacecraft are computed. The 
scene is then identified and the appropriate spectral factors are determin- 
ed for the final calculation of the radiances. Also, a bidirectional model 
R based on this scene type is used to account for anisotropy of the radia- 
tion in the estimate of the radiant exitance M for the pixel, i.e. the field 
of view of a single measurement, by use of equation .(1): 

M = it L/R (1) 

The scene identification procedure for ERBE first classifies the under- 
lying surface by its geographic type (geotype): land, ocean, snow, desert 
or land/ocean mix (coastal regions). The cloud cover is defined in terms 
of four classes: clear (less than 5% cloud cover), partly cloudy (between 
5% and 50% cloud), mostly cloudy (between 50% and 95% cloud) and overcast 
(more than 95% cloudy). 

The cloud class for a given measurement pair (shortwave and longwave 
radiances) is selected according to the measurement pair in a manner similar 
to a threshold method, but modified to account for the variance associated 
with each class, as illustrated in Figure 3. There is a distribution of 
.measurement pairs within each class, as represented by the light lines. The 
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shortwave-longwave plane is partitioned into regions according to which is 
the most probable scene identification for the measurement pair. The bound- 
aries of these regions are indicated by the heavy solid lines in Figure 3. 
Thus, if a measurement pair falls within region I, given that the surface 
is ocean, it is identified as "clear ocean". For nighttime scene identifi- 
cation, the cloud class assignment is made according to the longwave radi- 
ance only, in a manner analogous to the daytime case. 

B . Simulation approach 

It is seen that the ERBE scene identification algorithm is a bispectral 
technique using SW and LW measurements. In order to simulate the operation 
of the scene identification algorithm for a given scene, it is necessary to 
compute the SW and LW measurements for the scene. The measurement for a 
radiometer channel is given by 

co 

m c = / S, c L, dX (2) 

0 

where S^ c is the spectral response of channel c. In this paper, rad- 
iative transfer theory is used to compute spectral radiances for a variety 
of scenes, and the resulting measurements will be computed for the ERRE 
scanning radiometer channels and for the AVHRR channels using equation (2). 

Shortwave and longwave spectral radiances are computed for a number of 
model atmospheres and surfaces. In order to avoid excessive computational 
requirements, a two stream approximation is used to compute the shortwave 
spectral radiant exitances. The spectral radiances are then calculated by 
use of spectral bidirectional models R\ which were computed by Stuhlmann 
et al . (1985) and stored on tape. These SW and LW spectral radiances 
constitute the Level I data. The measurements were then computed for a 
given radiometer channel by use of these spectral radiances and its spectral 
response in equation (2), and make up the Level II data. Finally, the means. 
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variances and correlation coefficients are computed for this set of SW and 
LW measurements. These statistics are the Level III data set, which is re- 
quired for the scene identification algorithm. Because the model atmospheres 
and surfaces for the radiance calculations were not statistically distributed, 
this procedure also does not yield a statistically valid distribution of 
scenes and their variances, which would be desireable. However, it does 
provide reasonable numbers for the shortwave and longwave means, their var- 
iances and correlations, which are needed for the scene identification al- 
gorithm. The reasonableness of the approach is demonstrated by comparison 
of the Level III results with ERBE and AVHRR measurements. The AVHRR meas- 
urements are from the V5 tape, which is an ERBE produced magnetic tape for 
validation. This tape contains concurrent measurements from the ERBE and 
AVHRR instruments on board the NOAA 9 spacecraft. In order to provide a 
truth set for comparison with the ERBE scene identifications, the maximum 
likelihood method is adapted to the AVHRR data. In addition, the method 
of Eyre et al . (1984) is used with the AVHRR data for cloud discrimination 
to provide assurance that the AVHRR scene identifications are reliable. 

It is then determined which AVHRR pixels fit within a given ERBE pixel, 
and the probabilities are computed for the scene identification of the 
ERBE. Finally, the probability matrices for the scene identification are 
used for the evaluation of the errors which are incurred in data inversion 
(Smith et al . , 1984). The approach is shown schematically in figure 4. 
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III. COMPUTATION OF THEORETICAL FILTERED RADIANCES 

In the computation of Lx , it is assumed that radiation in the short- 
wave spectral range is reflected solar radiation and that radiation in the 
longwave is Earth-atmosphere emitted radiation only. 

Ill A. Radiative Transfer Computations 

For LW radiation, scattering is neglected and the surface is assumed to 
be perfectly absorbing. Thus, the spectral radiance along a ray at a given 
zenith angle is given by integrating the absorption and emission along the 
ray. Computations are made for each of nine zenith angles: 0°, in 0 ,...,8D°. 
Because solar radiation does not contribute in this spectral range, there 
is no dependence on 0g and <)>. Spectral radiances are computed as 
described by Schmetz and Raschke (1981). Briefly, the LW model computes 
spectral radiances in each of 51 spectral intervals, ranging from 4 to 
400 urn. Absorption and emission by H2O, O3, H2O polymers and "uniformly 
mixed gases" (i.e. CO2, CH4, N2O) and in clouds, water droplets are taken 
into account. Longwave radiant exitance is next computed from these spectral 
radiances by integrating over the zenith angle and the spectral ranges. 

For the computation of SW radiances, in the range 0.20 to 3.7 u m, mul- 
tiple scattering must be considered. In order to reduce the computational 
requirements, spectral radiant exitances are computed by use of the two- 
stream approximation and the spectral bidirectional models of Stuhlmann et 
al . (1985) is applied to all SW wavelengths to get spectral radiances. 

The computation of SW spectral fluxes includes Rayleigh and aerosol scatter- 
ing and absorption by water vapor, aerosols, ozone, the "uniformly mixed 
gases" (i.e. CO2, CH4, N2O) , and in clouds, water droplets. Transmission 
functions for the atmospheric gases are taken from Kerschgens et al . (1976), 
and aerosol radiative properties are those recommended by IAMAP (McClatchey 
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et al . , 1977). The spectral reflectivity of the lower boundary was treated 
as by Kerschgens et al . (1976). 

Cloud was considered to be a one kilometer thick homogeneous layer 
within the atmosphere, with variable cloud top heights from one to 15 km. 

The cloud droplet size distribution is the Sc-2 from Schmetz et al . (1981) 
as measured in stratocumulus clouds, with an effective droplet radius of 
7.9 ym. The radiative properties, including single scattering albedo and 
asymmetry factor are the same as those computed by Stuhlmann et al . (1985), 
using Mie theory. The optical depth at X = .55 ym was determined from the 
droplet size distribution, and the optical depth at other wavelengths was 
then computed, based on Schmetz et al (1981). Varying the cloud optical 
depth between 1 and 100 results in a range of column liquid water content 
of 5-500 g-m"2. 

In order to compute shortwave spectral radiances from the spectral rad- 
iant exitances which were based on the two-stream approximation, spectral 
bidirectional models of four surface scenes and two clouds (h c i = 3 km and 
9 km.) by Stuhlmann et al . (1985) were used, which are based on an iterative 
scheme by Raschke- (1971). The calculations by the 2-stream approximations 
used in the present study are based on the same atmospheric input data. 

The computations by Stuhlmann et al . (1985) treat the lower boundary as 
Lambertian except for water surface, for which the reflection is considered 
to be partly diffuse and partly specular. The specular part is computed 
from Fresnel's formula, and the surface is assumed to be wavy. 

Ill B. Scene Description 

Atmospheric profiles of temperature and humidity were taken from 
McClatchey et al . (1977) for five climatological conditions: subarctic 
summer SAS, midlatitude winter MLW, midlatitude summer MLS, subtropic summer 



- 9 - 


STR, and tropics TRO. These profiles are shown in figure 5. Table I shows 
as an example the content of the data set of one atmospheric set of one 
profile (subtropic summer) for both the SW and the LW range. The aerosol 
density is described by the visibility parameter 3. The aerosol profile 
of McClatchey et al . (1972) is used; for this case, 3 is 56 km. A value of 
3 = 23 km was also used, corresponding to an increased aerosol density. 

In addition to the surface temperature of the model, surface temperatures 
of T m0 del“ c and T model +c were used 1n the radiance computations. This pro- 
vides a simulation of the effect of diurnal cycle on the radiances. The 
change c was taken to be 4°K for SAS, MLS and MLW and to be 6°K for STR and 
TRO. All surface scenes are sorted into four classes according to the ERBE 
classifications: ocean, land, ice/snow, desert. 

Spectral radiances were computed for each angular bin for each model 
type listed in Table II. Thus, for each surface type there is a range of 
cloud heights and optical thicknesses, in addition to. the three surface 
temperatures. The set of spectral radiances thus computed constitutes the 
"level I" data in figure 4. Given the spectral response of various satellite 
instruments, the resulting measurements of these scenes can be computed by 
use of equation (2) . 

Ill C. Spectral Responses of Spacecraft Instruments 

The spectral response of a satellite radiometer to incoming radiation 
is determined by its optics and the nature of the sensing element. 'For the 
total channel of the ERRE scanning radiometer, a Cassegrain telescope con- 
sisting of two front surface aluminum mirrors focuses the radiation onto a 
bolometer, which measures the total radiation impinging on it. The spectral 
response of this channel is thus the product of the spectral reflectivity 
of each of the mirrors and the spectral absorbtivity of the coating of the 



- 10 - 


bolometer, and was shown in figure 1. The shortwave channel is the same, 
except for a suprasil filter in front of the detector, which effectively 
blocks longwave radiation but passes shortwave radiation. The longwave 
channel has a diamond filter in front of the detector, which blocks short- 
wave radiation but passes longwave radiation. The spectral response of 
each of these two channels is that of the total channel times the spectral 
transmission of the filter, and was shown also in figure 1. 

Because data from the ERB scanning radiometer on the Nimbus 6 and 7 
spacecraft and the Advanced Very High Resolution Radiometer AVHRR aboard 
the NOAA 9 spacecraft are used, it is necessary to consider the spectral 
responses of these instruments. The shortwave and longwave channels of 
the ERB scanning radiometer on the Nimbus 6 and 7 spacecraft have practic- 
ally the same spectral responses as those of ERBE. Data from channels 1, 2 
and 4 of the AVHRR are also used, and their spectral responses are shown 
in figure 6. 

The measurement for a given radiometer channel was thus computed by use 
of its spectral response and the spectral radiance of the scene in equation 
(2). These simulated measurements constitute the level II data for this 
study. 

Ill D. Means and Variances of Filtered Radiances 

For each scene type, spectral radiances and simulated measurements were 
computed for each of three surface temperatures and a range of cloud top 
heights and optical thicknesses for each viewing zenith, solar zenith and 
azimuth angle, as shown in Table II. For each scene type and angular bin, 
average radiances, standard deviations and correlation coefficients between 
the shortwave and longwave measurements were computed. These data, which 
are called the level III data set, are the kind which are required for in- 
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put to the scene identification algorithm for ERRE data processing. 

Typical results for the ERBE scanning radiometer are shown in figures 
7a and 7b, respectively. Average shortwave and longwave channel simulated 
measurements and their standard deviations are shown versus viewing angle 
bins in the backscatter direction (150° < (|> < 180°) and solar zenith angle 
in the range 32.5 to 37.5°. The scene is tropical ocean for clear and with 
high and medium height clouds. The clear ocean scene shows the expected 
limb darkening in the longwave measurement and limb brightening in the short- 
wave measurement. In the longwave case the radiances of the high cold cloud 
(HK) increase with increasing 0 V because the atmospheric effects superim- 
pose the low cloud temperatures. For the same reason the shortwave radiances 
decrease at high viewing zenith bins. The simulation of the sun glint is 
another feature of the applied bidirectional models. (0° < <j» < 9°). 

The atmospheric temperature and humidity profiles (Me Clatchey et aT., 
1977) which were used for the longwave radiance computation do not consider 
the surface types which are used by the ERBE scene identification algorithm. 
For the present study, the effect of surface type on the longwave measure- 
ment is accounted for by empirically determined factors, which are obtained 
from mean climatological LW exitance data. Thus, the longwave radiance for 
an atmospheric profile p and scene ss is computed by 

Lp LW ( 9 v’ ss ^ = L p LW ( 9 v) x V ss ) (3) 

where _ 

fp(ss) = M*(tp,Xp)' (4) 

2 n f T lw (9 ) sin 9 cos 9 d 9 
0 

ss = surface scene (ocean, land, ice/snow, desert) 
p = atmospheric profile (subarctic winter,..., tropic) 

TC*(tp,x p ) is a mean value of all latitude bands (xp) and the season (tp), 
which corresponds with one of the five available atmospheric model profiles. 
The factor fp(ss) is mainly in the range .9 to 1.1. 
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IV. APPLICATION OF THEORY TO ERBE SCENES. 

Simulated measurements were computed for the ERBE SW and LW scanner 
channels using the spectral responses for these channels together with the 
computed spectral radiances. The means, variances and correlation coeffic- 
ients were computed for these simulated measurements. The 1-sigma ellipses 
are shown in the LW-SW plane in Figure 8. The locations of the centers of 
the ellipses in this plane are given by the means and their sizes and 
orientations are determined by the correlation coefficients and variances. 
Also shown in this figure are real ERBE measurements, with plotting symbols 
indicating scene types as determined by the ERBE scene identification al- 
gorithm. The area is located at the Pacific west of California and Central 
America, for April 4, 1985. Here only "clear ocean" and "overcast" scenes 
were chosen. 

The ERBE determined overcast scenes are in good agreement with the 
model "high clouds", with a cluster near the clouds with medium optical 
thickness. This result is also confirmed by AVHRR data, as will be discuss- 
ed later (see VB). On the other side the "clear-ocean" pixel cluster is 
"cooler" and a little "brighter" than the pixel average. This might be due 
to the large ERBE pixel size (44 km *65 km at nadir), which most probably 
leads to cloud contamination. In fact an ERBE scene is denoted as "clear" 
when the pixels contain less than 5% cloud, which of course increases the 
SW and decreases the LW radiance. 

Another obvious reason for possible disagreements is the necessary 
limitation in model input data (albedo models, visibilities, surface temp- 
eratures, cloud parameters, atmospheric models). This modifies the mean 
radiances and diminishes their variances. In spite of this, the data set 
provides agreement which is good enough that the present method is seen to 



- 13 - 


provide an acceptable tool for scene identification. 

One conjectures that the data points between the "clear" and "overcast" 
may be due to high thin and high medium thick clouds. However, they may 
also be due to partly filled fields of view, or partly cloudy scenes. Use 
of higher resolution data such as AVHRR imagery is necessary to distinguish 
between these alternatives. 
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V. ERRE-AVHRR (NOAA-9) SCENE-ID INTERCOMPARISON 
V A. Technique 

The Advanced Very High Resolution Radiometer (AVHRR) on board the NOAA-9 
spacecraft with the ERBE instruments provides an excellent opportunity for 
an independent check on the scene identification as computed by the ERBE 
procedures. The smaller field of view of the AVHRR gives more pixels which 
are clear or cloud filled than does the ERBE scanner. Also, techniques have 
been developed for the identification of clouds with AVHRR data. In this 
study, the maximum likelihood estimator (MLE), which is the ERBE scene iden- 
tification algorithm as discussed earlier, is applied to the AVHRR data to 
compute scene types to provide a comparison with the ERRE computed scene 
types. Next, the AVHRR scene types are computed by use of the method of 
Eyre et al . (1984) and comparisons are made. 

Data for these intercomparisons are taken from ERRE V-5 data tapes, 
which contain data from both the ERBE and AVHRR radiometers for selected 
regions of the Earth, for April 14, 1985. Two test sites were chosen, as 
indicated on the map in figure 9. The first is located in the Pacific 
Ocean, west of and Central America, as the spacecraft moved from 4.9°N, 

132. 4°W to 9.7°N, 133. 6°W, covering approximately region number 1197 of 
the ERBE 2.5°x2.5° grid as the scanning radiometer operates in across-track 
mode. The second is in the Atlantic Ocean west of Africa but including 
some coastal areas, as the spacecraft moved from 4.9°N, 4.7°W to Q.7°N, 
5.9°W, covering approximately region number 1219 of the ERBE grid. 

The MLE method selects a cloud class for a given pair of measurements. 

In doing so, it treats each pixel independently to attach one of the a 
priori modeled scene types, characterized by radiance averages, standard 
deviations and correlation coefficient. In order to apply this method to 
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AVHRR data, these statistics, which we call level III data, must be avail- 
able for the channels of the AVHRR, The methods described in sections II 
and III were used to generate this required data base, for the 20 scene 
types. 

The present study uses the simplifying assumption that the scene is 
either clear or completely cloudy. Arking et al . (1983) mention the im- 
pacts of this error and use a "maximal clustering technique" to determine 
cloud fractions of the field of view of each pixel. A different method 
more suitable for this study is described later (§VD). Another feature of 
the MLE method is a quadratic dependence on the standard deviation. Pixels 
which are much warmer than the modeled clear scene LW value with its norm- 
ally rather small variance (for example, clear "desert" or "land") could 
thus be identified as a cloud (significantly higher variance). This "wrap 
around" is avoided in this study by including some restraints in the MLE- 
probability computation, as is done in the ERBE scene- identification algor- 
ithm. 

In order to apply the MLE method for scene identification from AVHRR 
data, a new level II data set suitable for the AVHRR is computed using 
AVHRR channel 2 and 4 filter functions (NOAA Users' Guide), centered around 
.9 and 11. urn, respectively, to simulate SW and LW measurement pairs. A 
level III data set for the AVHRR measurements was created which includes 
average filtered radiances for the same 2f> scenes listed in Table II. 
Channel 2 was chosen instead of channel 1 (.55-. 75 ym) because in this 
channel calibrated radiances of various scenes cover a broader range than 
the respective radiances of channel 1, e.g., clouds appear brighter and 
clear scenes appear darker in channel 2 than in channel 1. This behavior 
obviously allows a better scene discrimination in the SW-range. 



- 16 - 


The counts of channels 2 and 4 were calibrated and transformed into rad- 
iances using the calibration coefficients stored on the V-5 tape using the 
procedure described by Rossow et al . (1985). Percent scaled radiances for 
each channel are computed from the counts using tabulated gains and offsets. 
The effective solar irradiance for each channel is calculated by integrating 
the product of the spectral response function for the channel and the solar 
irradiance, as tabulated by Neckel and Lab (1984). The radiance is then 
the percent scaled irradiance multiplied by the effective solar irradiance 
for each channel . 

With the geographic information of each pixel, an underlying surface 
scene is taken from an external data set. The a priori model values of 
this clear scene and of a variety of model clouds above this surface together 
with the original AVHRR measurement pair serve as input data in the MLE 
procedure. The model scene with the highest probability is regarded as 
the "true" scene for the AVHRR pixel. Uncertainties -in this scene deter- 
mination are due to the presupposition of a Gaussian distribution of the 
pixel in the SW-LW plane and of course to the fact that the model scene 
possibly does not fit the reality. 

One-sigma ellipses in the LW-SW plane for the level III data set for 
AVHRR measurements are shown in figure 10 a and b for tropical ocean, land 
and clouds. Calibrated LW and SW AVHRR data (channels 2 and 4), indicated 
as dots, are compared with'the model scene values as in figure 8. These 
data are AVHRR GAC measurements from the NOAA 9 satellite on 14 April, 1985 
for the two regions earlier described. The GAC data from the V-5 tapes are 
much smaller (~12 km^) than the ERBE pixel (44x65 km, N0AA9). It is 
seen that the AVHRR data lie in the envelope of 1-sigma ellipses of radiances 
as computed by radiative transfer. While the first case a) consists of cold 
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ocean scenes with various, cloud types, the second test area contains clear 
land pixels and only a few (warm) ocean pixels. The thin cloud types over 
land are significantly brighter than the respective clouds over ocean in a) 
and are in good agreement with the modeled scenes (e.g. scene types 9,9 and 
10 in Figure 10 b.). The influence of the surface can be neglected for med- 
ium thick (scenes 13,14 and 15) and thick (scenes 19,20) clouds. This agree- 
ment lends support to acceptance of the identification of these scene types, 
although as before it is noted that it is assumed that all fields of view 
are either full or empty. 

The following problems occur. The variance of the SW radiances for the 
ocean is small compared to that for partly cloudy over ocean. As a conse- 
quence, the MLE method will tend to identify many scenes as cloud. Also, 
it is difficult to discriminate between clear ocean or clear land and thin 
and low clouds, as their radiative effects are small. Fortunately, for this 
reason the influence of these scene identification errors is small. Finally, 
many GAC pixels which are classed as clear may well be cloud contaminated. 

V B. Spatial Colocation of ERBE and AVHRR pixels 

In order to compare the scene identification results from the AVHRR with 
those from the ERBE scanner, it is necessary to colocate them spatially. 

The ERBE scanning radiometers have an optical mask in front of each detector 
which forms a hexagonal instantaneous field of view with the 3° width in the 
scan direction and 4.5° width in the perpendicular direction. The shape of 
the ERBE pixel is also modified by the time response of the instrument as it 
scans. Because of this the shape of the surface area covered by an ERBE 
pixel is approximated in this study hy superimposing an ellipse and a rect- 
angle with the footprint of the pixel in its center, as shown in figure 11. 
The length of the major axis "a" along the scan line is approximated as half 



- 18 - 


of the distance of two neighboring pixel -footprints. The minor axis "b" is 
half of the distance of two neighboring footprints along the nadir and is 
held constant along the scan line. With the geometric information of both 
the ERBE and the AVHRR data the distance of the AVHRR footprint perpendicular 
to the ERBE scan line "d" and the distance to the edge of the ERBE pixel is 
determined and compared (figure 11). If |d|<b', the AVHRR pixel is regarded 
as being within the ERRE pixel. 

V C. Method of Cloudiness Oetermi nation 

As mentioned above, it is necessary to account for partially cloud-filled 
AVHRR pixels in order to estimate an average cloudiness of the total area. 
This was done by applying a technique which makes use of the different 
behaviour of AVHRR channels 1 and 2 reflectances (ri.r?). For ocean scenes 
ri is larger than r2 and for land and desert the opposite is the case, which 
might be due to differing aerosol scattering in the channel 2 range or 
to the spectral characteristics of the incoming rays.. For the high resolu- 
tion AVHRR-LAC data (1 km^), Kriebel (1986) and Saunders (1986) found that 
for clear ocean r2/rl ~ 0.5 and for land r2/rl > 1. The maxima for clouds 
were centered at 0.9. This value can vary between 0*7 and 1.0 according to 
cloud top altitudes as a result of change of ozone and water vapor content 
above the clouds. 

In order to apply this behaviour to identify partly cloudy AVHRR pixels, 
as a first step one of the 20 pre-defined scene types is attached to each 
AVHRR pixel with a Maximum Likelihood Estimation method. Also the ratios of 
the reflectances ^/r^ are determined for each pixel and the corresponding 
scenes are sorted into histograms as shown in figure 12. Because of the 
larger pixel sizes the peaks are not so well determined, which points to a 
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probable cloud contamination of the clear scenes. Clear ocean, desert and 
land scenes are clearly separated from each other and show overlapping areas 
with the cloud cluster. All pixels with ratios in these areas where 

regarded as cloud contaminated or as not totally cloudy in the case of pre- 
defined cloud scenes. Boundaries were chosen for clear, mixed and cloudy 
scenes depending on the respective condition of the histograms of each test 
area. Then the AVHRR pixel cloudiness amount c is determined. It is 0 for 
clear scenes and 1 for totally cloudy pixels and linear interpolated between 
0 and 1 if ^/r^ is beyond the predefined boundaries (e.g. for an ocean scene 
scene in figure 12, 0 < c < 1 if 0.75 < < 0.9). 

As a final step these cloud amounts c of all AVHRR pixels belonging to 
one ERBE pixel results in a mean percentage cloudiness ~c of the pixel area 
which is according to the ERBE scheme : 

0 % < c < 5 % : clear , j = 1 

5 % < c^ < 50 % : partly cloudy. , j = 2 

50 % < "c < 95 % : mostly cloudy , j = 3 

95 % < "c < 100% : overcast , j = 4 

The type of cloudiness is thus computed for each area covered by an ERBE 
pixel (see Table II). These scenes j, here regarded as "truth", can finally 
be compared with the respective ERBE scenes of the same area. 

The ERBE-AVHRR scene intercomparison were done for all pixels of a test 
area and the results were sorted into matrices as shown in Table III for 
region 1197. In part (a) of Table III, all intercomparisons are taken into 
account. In 353 cases the ERBE scene is computed to be clear ocean (i=l) 
whereas partly cloudy scenes are determined from the AVHRR measurements 
using the method described above. The example shows clearly some other fea- 
tures that are found. In the clear and partly cloudy cases for AVHRR, there 
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are deviations of the ERBE to both more and less cloudiness. A slight tend- 
ency for underestimation can be seen at partly cloudy scenes. Yet this 
feature is very clear at mostly cloudy and overcast scenes, where ERBE for 
instance determines 687 mostly cloudy scenes that were classed by the AVHRR 
procedure as overcast. 

The matrices in part b) of Table III contain all those intercomparisons 
with the restriction that a single AVHRR determined cloud type is dominant 
for each cloud type (c (cloud type) > 50 % ). From the first matrix it is 
evident that the underestimation is due to thin clouds and in a less degree 
to low and middle high clouds. On the other hand, there is almost always a 
fairly good agreement at medium, thick and high clouds. In all cases only 
very few overestimations of cloudiness from the ERBE are observed. 
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VO. Figures of Merit for Intercomparisons 

Two methods are chosen as means to compare the matrices of different 
areas, from now on called probability-matrices P^j according to Smith et 
al . (1984). As the first characteristic the sum of the diagonal elements 
(denoted -3,. ..,3) as shown below are determined and their percentage of 
the matrix-elements total is plotted as histograms (Figure 13) for the five 
cases which were studied. 
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In order to characterize the degree of misinterpretation of scenes a weight- 
ing matrix is defined as 


( a ij)i*l,4 = 
j-1,4 


12 3 4 

2 12 3 

3 2 12 

4 3 2 1 


Two more matrices are necessary in order to define a parameter which is a 
measure for the distribution of the elements within a matrix : 


( c ij) = { 


a ij 


Pij > 0 
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n , Pij * o 

Also the following quantities are defined: 
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With these numbers and denoting the number of rows/columns as n, the required 
parameter MM can be defined : 


T = 
S = 


1 A 
n-T 7 
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7TT 
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1 < A < n , 0 < T < 1 

P 

l<C_<n , 0 < S < 1 

0 


n- 

Def. : MM = T /T 

MM stands for "Matrix Moment" and is 1 if all elements are on the middle 
diagonal (absolut agreement) and 0 if they are located in the lower left and 
upper right corners (total disagreement). T(A/P) is the normalized average 
weighting factor of the matrix Pjj with the worst result at A/P = 4. C/D 
is a factor that takes the distribution of the elements within the matrix 
into account. Figure 13 demonstrates that for equal values of A/P, MM de- 
creases if the matrix-elements are changed toward the corners. 

The following example shows two matrices with 14 elements in different 
positions but with the same ratio A/P = 24/14. 
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MM = 0.574 


Although there are more elements on the middle - diagonal at M2 MM is 
lower because of three elements being at the worst position in the lower 
left corner. 
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V E. Discussion of Intercomparison Results 

Figure 14 shows histograms of the normalized diagonal sums and the 
matrix moment MM for several cases. The histograms show similar patterns: 

MM between is between .65 and .73 in all cases and there are only few strong 
deviations from the middle diagonal. In four cases there is generally more 
underestimaion of cloudiness according to the ERBE scene identification 
method. This underestimation occurred often at thin and low clouds, whereas 
there is almost always a good agreement at medium, thick and high clouds. 
Thus ERRE tends to interpret thin clouds in lower levels as clear scenes 
(at high cloud amounts), or more accurately, the AVHRR scene identification 
method determines more thin and low clouds in spite of its tendency to de- 
note those pixels as mixed scenes (see VD, 0 < c m -j x < 1). 

A significant overestimation of cloudiness is found in only one case 
(April 14, 1985 for 1219 Atlantic/Northwest Africa). This is the only test 
area that contains about 30% clear and cloudy land, desert and coastal 
scenes according to the a priori ERBE surface scene map. Reasons for this 
might be the fact that the "land" scenes of this part of Africa are found 
to be very hot and sometimes as bright as "desert" scenes. Together with 
the "wrap around" effect of the MLE method this could be one reason for 
about 20% of all clear land scenes (AVHRR) being denoted as "partly cloudy" 
by the ERBE scene method. Finally there are several bright and extremely 
hot desert scenes determined by the AVHRR in areas the ERBE geotype map 
has predefined as "coastal" areas. Due to the "wrap around" effect these 
obviously "clear" scenes are misidentified as partly cloudy and in several 
cases even as overcast scenes by ERRE. 

On April 19, 1985 at the same test area (1219) there are several wrong 
(overcast) AVHRR scenes found in the sun-glint area of clear ocean where 
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the ERBE quotes the scene rightly (according to albedo and surface tempera- 
ture) as cloud-free ocean. That indicates a reasonable handling of the 
sun-glint in the ERBE inversion algorithm. This statement is also confirmed 
by the fact that neglecting all comparisons in the sun-glint area does not 
improve the results of the probability matrices. 

The matrices shown as an example in Table III were also sorted into all 
combinations of angular bins. It can be seen that misinterpretations of 
ERBE scenes are at no preferred angular bin, e.g. there are no better inter- 
comparison results in backward than in forward scatter direction. Finally, 
if the ERBE determines a scene as overcast, it is most probably true because 
there are only very few overestimates at high cloud amounts (overcast 
instead of mostly cloudy) for all cases. 
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VI. ERROR PROPAGATION IN THE ERRE INVERSION ALGORITHM 
VI A. Theoretical Considerations 

In this study, the approach of Smith et al . (198A) will be used. The 
total error of the radiant exitance for each pixel is the sum of two basic 
components. The error in scene identification results in discrete changes 
in the spectral correction and bidirectional functions (Fig. 15), or branch 
ing in^ the computational flow. The second error component contains term 
errors (bidirectional function- and spectral factor- variations, instrument 
noise) and can be treated in a linear manner. In this study some first 
estimates for the errors due to scene misidentification will be presented. 

Let a scene which is in truth type "j" be identified as type "i". Then 
the expected error in the estimate of radiant exitances M is defined as 


E( e Mi ) * Z e Mij Pi j (5) 

J 

where P-j j is the probability that a scene is in truth type j and that it 
is identified as type i. This error, if nonzero, represents a bias in the 
analysis of the data. Assuming that this bias is zero the variance of 
this error is 

V ^ e Mij) = ? e Mi j P ij ^ ^ e Mi j P ij ^ = f e Mij p ij ^ 

The probability Pij is the product of the conditional probabilities that 

a scene is in truth type j and that it is identified as type i. As such, 

it can be written as 


p ij = p j / n j( x ) dx (7) 

R i . 

where x is a two dimensional radiance space (L SW ,L IW ). The probability of 
occurence of various scene classes Pj represents the statistical makeup of 
the cloud field, and Hj(x) is the probability that if the scene is of type 
j that it will have value x. The domain of integration is the part of the 
x or (L SW ,L^ W ) plane in which the scene identification algorithm identifies 
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the scene to he of type i. The probability distribution H j ( x ) will be 
considered to be a normalized probability distribution of the measured 
radiance pair (L sw ,lJ w ) for scene type j and can be approximated as 

e-° (A) 


V x> * < °sw °lv. 2 ’ 




0 ■ ( t * 4 


- 2 p 


u v 


°sw a lw 


) 


u = LSW _ L sw 
v = lJ w _ iJ w 

The means, variances and correlation coefficients of shortwave and longwave 
radiances are considered to be known for each scene type j. The integration 
is over the region which is defined by the scene identification algorithm 
using equation (8) as being scene type i, indicated in figure 16 by the 
cross-hatched area. These scene areas are delineated by ellipses which are 
given by the equation 


Pi Hi- Pj Hj (9) 

Because the Gauss distribution for scene classes may be an unrealistic 
assumption for points well away from the centroid of each scene type, it 
is desirable to use the Pij as computed in the last chapter. The AVHRR 
scene type is again considered to he truth j. 


AVHRR/j 

ERBE/i 

1 

2 

3 

4 

1 

n ll 

"12 

"13 

"14 

2 

n 21 

"22 

"23 

"24 

3 

"31 

"32 

"33 

"34 

4 

"41 

"42 

"43 

"44 

Nj : 

Ni 

N 2 

n 3 

n 4 


= N tot 



- 27 


The portion of scene class j identified as i (the cross-hatched area in 
figure 16) is thus evaluated empirically rather than by the integral ex- 
pression in equation (7): 


Ni 


\j "1j n i j 

x = 

Ntnt TH 


ij 


(10) 


'•tot "j 'Hot 

The ratio of Nj and the total number of pixels in the test area Ntot an 
empirical evaluation of P j , and the ratio of n-jj to Nj is an empirical 
evaluation of P-jj. The P-j j are only valid for a certain area and based on 
a "truth" which of course is an estimate, based on many uncertainties, 
by itself. 

If the scene is identified as type i the shortwave and longwave radiant 
exitances at TOA, expressed as components of a two-dimensional vector, are 
estimated as 


Mi ■ it Qi Cj m (11) 

whereas the correct exitance for scene class j would be 

Mj = u Qj Cj m (12) 

Q is a 2 x2 diagonal matrix with 1/R on the diagonal for the shortwave- 
and longwave component, respectively. C is the 2x3 spectral correction 
matrix and m a vector of ERBE shortwave, longwave and total channel radiances. 

The error eM-jj in the estimate of radiant exitance due to scene 
errors is then 


£Mij = 'w (Qi - Qj)Ci!Il + ifQj(Ci - Cj)iD 
The covariance matrix fore^jj is given by (analogous to (6)): 


( 13 ) 
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cov(e^-j (scene)) = 

* 2 Z (Oi - Oj) lyij + cov^)] (0| - oj) 

J ^ 


if 2 I P^j Oj (C^ - Cj) + cov(mj)] (cj - cj) oj 

J ^ 


+ ^ j.jij ^ K i j + K ij^ 

Here the approximation is made that 


( 14 ) 


P • J n • (x) dx m m T 
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p ij ^Eli Ell + cov(mj)) 


(15) 


The second moment of the distribution is in this way evaluated by use of 
the a priori mean values of the identified scene radiances. 

The first term on the right hand side of equation (14) is the error 
variance due to bidirectional modeling error, the second term is the error 
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cation error, and the remaining parts are coupling terms between these two 
effects. 


Preliminary Estimates of Exitance Errors 
The shortwave and longwave exitance standard deviations are the square 
roots of the diagonal elements of the covariance matrices. For each pixel 


the percentage deviation can be calculated: 

a i 

p, = — • 100 [%] (16) 

Mi 

Figures 17 and 18 show preliminary results of these errors for the ERBE 


scenes in region 1197 on April 4, 1985. The probability matrix for this 
case is: 
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Because of the low likelihood of underestimating an overcast scene the 
shortwave errors for this scene type form the minimum areas in figure 17. 
Error maxima are almost always connected with mostly cloudy and sometimes 
partly cloudy scenes. Reasons for this are the great amount of variability 
in cloud types, altitudes and total amount (5-50 and 50 - Q5 %) which 
causes the radiometer signal m and finally the radiance estimate M_. The 
extremely high deviations, indicated with "H" in figure 17 are only very 
scattered and mostly caused by overestimation of cloudiness when compared 
with corresponding "true" scenes. Table IV summarizes these results for 
one case. Errors are averaged over three parts of an ERBE scan line along 
the flight track. Although these values are based on a limited statistical 
sample and only valid for a certain area, they demonstrate some typical 
features. Obviously the longwave errors are negligible for all scene 
classes and also the shortwave deviations are very low for overcast scenes. 
For all scenes there is an increasing absolute error. in forward scatter 
direction but only for clear ocean this is also the case for the percentage 
error. For all cloudy scene types the shortwave exitances are more increas 
ing than the absolute errors along the scan-line toward the forward scatter 
ing of the sun light, and therefore the percentage maxima are in the nadir 
region. This clearly indicates that the shortwave radiant exitance error 
due to scene identification error depends on the direction of the reflected 
sunlight for a clear ocean surface more than for cloudy scenes, because of 
the sun-glint effect over the ocean. 
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VII. CONCLUSIONS 

The scene identification probability matrices quoting the probability 
that a scene which is in truth type j is identified as type i are required 
for the evaluation of data analysis errors. This paper presents empirical 
evaluations of these matrices by checking the accuracy of the ERRE scene 
identification algorithm. This procedure attaches a scene (clear surfaces, 
partly and mostly cloudy above these surfaces and ocercast) to a shortwave- 
longwave measurement pair with a maximum likelihood estimation method. In 
order to perform an intercomparison of these scenes with an independent set 
of satellite data, a scene identification procedure is developed for AVHRR 
SW and LW data by using two radiative transfer models. These AVHRR RAC 
pixels are then attached to corresponding ERBE pixels and the results are 
sorted into scene or probability matrices. 

These scene-type intercomparisons show that there generally is a higher 
tendency for underestimation of cloudiness over ocean'at high cloud amounts, 
e.g. mostly cloudy instead of overcast and partly cloudy instead of mostly 
cloudy, of the ERRE scene identification relative to the AVHRR results. 

This is mostly caused by thin or low clouds, which on the other hand has 
only a small impact on the results because of the similarity of the respect- 
ive bidirectional functions. Another error source is the so-called "wrap 
around" effect in the maximum likelihood method which causes cloudiness 
overestimation over hot land and coastal scenes. 

Preliminary estimates of the errors in radiant exitances which are due 
to scene misidentification demonstrates the high dependency on the probabil- 
ity matrices. While the longwave error can generally be neglected the 
shortwave deviations have reached maximum values of more than 12 % of the 


respective exitances. 
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VIII. FUTURE RESEARCH REQUIRED 

It is planned as a continuation of this study to develop a nightime 
scene-identification method for AVHRR data using only the longwave channels 
3 (centered at 3.7 y m), 4 (centered at 11 y m) and 5 (around 12- m). The 
difference in brightness temperatures of channel 3 and 4 for example enables 
us to detect low clouds, fog or thin cirrus. This feature is described and 
used for image processing by Eyre et al . (1984) or for retrieval of cloud- 
free pixels at night by Saunders (1986). 

In order to estimate the total pixel error, the uncertainties in the 
signal (m sw , m-| w ), the filter factors (C) and the bidirectional functions 
(R) must be taken into account. These error sources propagate in the inver- 
sion algorithm in a linear manner. Further theoretical studies are necess- 
ary to estimate the variances of R sw and R] w . Also, the analysis of Smith 
et al . (1984) should be extended to the case of the expected error of a 
specific measurement. 
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8. Bispectral distribution of ERBE measurements aboard NOAA 9 spacecraft 

and modeled one-sigma ellipses within for a sample case. 

For April 4, 1985, area 1197 (see text). 

9. Location of orbit tracks for data periods. 

10. Bispectral distribution of AVHRR measurements aboard NOAA 9 spacecraft 

and modeled one-sigma ellipses within for a sample case. 

a) for April 4, 1985, area 1197 

b) for April 14, 1985, area 1219 (see text). 

11. Geometry of mapping AVHRR pixels into ERBE pixels. 

12. Histogram of ratio Channel 2/Channel 1 for AVHRR data for differing 

scene types. 

13. Dependence of matrix moment MM on A/P and C/D. 

14. Histograms of sums of matrix diagonals 

15. Branching of computational flow due to scene identification. 

16. Schematic of scene identification probability computation. 

17. Standard deviation of shortwave radiant exitance errors due to scene 

identification errors. 
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TABLE I 

Example of one level I data set 
Model: Subtropic - summer 

a) Range: 0.2 - 3.28' m (Shortwave) 

Spectral fluxes of each scene for all combinations of 

- 37 spectral intervals 

- 16 sun zenith angles (0,. . . ,72.5°) 

- 15 cloud top heights (1,...,15 km) 

- 12 optical depths ( 1 , . . . ,100) 

- cloud-free case 
Surface scenes 


Albedo model 


visibility (km) 

scene 

classification 

ocean 


56 

ocean 

ocean 


23 

ocean 

limestone. 

light soil 

56 

desert 

dry meadow 

, grass 

56 . 

land 

savannah 


56 

land 

sand 


56 

desert 

sand 


23 

desert 

wasteland. 

subtropic low 

56 

land 

wasteland. 

subtropic middle 

56 

desert 

wasteland. 

subtropic high 

56 

desert 

wasteland. 

subtropic high 

23 

desert 


b) 4.09 - 400 y m (longwave) 

Spectral fluxes for all combinations of 

- 51 spectral intervals 

- 3 surface temperatures (T m0 deT c » T model ♦ T model +c ) 

9 viewing zenith angles (0,...,80°) 

- 15 cloud top heights (1,...,15 km) 

- 9 optical depths (1, 2, 4, 6, 8, 10, 12.5, 15, 20) at y = .55 y m. 

- cloud-free case 
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TABLE II. 

MODEL SCENES 

Model Surface Cloud Cloud 


Number 

1 

Ocean 

Thickness 
cloud free 

Hei ghts 


2 

Land 

cloud free 

- 


3 

Ice/snow 

cloud free 

- 


4 

Desert 

cloud free 

- 


5-7 

Ocean 

thin 

low, middle, 

high 

8-10 

Land 

thin 

low, middle. 

high 

11-12 

Ice/snow 

thin 

middle, high 


13-15 

Ocean/land/desert 

medium 

low, middle. 

high 

16-18 

Ice/snow 

medium 

low, middle. 

high 

19-20 


thick 

middle, high 



Model numbers 5-12: low 1 < h < 3 

middle 4 < h < 8 

high 9 < h < 15 

13-18: low 2 < h < 4 

middle 5 < h < 9 

high 10 < h < 15 

19-20: middle 5 < h < 10 

high 11 < h < 15 


Cloud thickness 

Optical thickness 

thin 

1 < T < 6 

medium 

8 < t < 20 

thick 

30 < t < 100 

Angular bins, degrees 

: 

0q: 0-2.5, 2.5' 

-7.5, 7.5-12.5, ..., 67.5-72.5. 

CD 

< 

0 

1 

H— * 

o 

• 

• 

• 

V* 

• 

G 

cc 

1 

o 

<|> : 0-9, 9-30, 

30-60, 150-180. 
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Table III 

Scene matrices for : 
region 1197 April 14, 1985 

AVHRR : j=l,2,3,4 scene type from averaged cloudiness c 

ERBE : i=l,2,3,4 " " " ERRE scene i.d. 

a) all pixel intercomparison : 


AVHRR/j 

1 

2 

3 

4 

ERBE/i 

1 

543 

353 

74 

0 

2 

177 

688 

390 

25 

3 

32 

283 

803 

687 

4 

0 

1 

45 

811 


b) 

more 

than 

50 

% of all 

AVHRR 

scenes covering 

one ERBE 




pixel 

(mostly cloudy. 

overcast) 

are : 


- 




1. thir 

i clouds 



2. 

medium clouds 

3. 

thick 

clouds 

AVHRR/j 

ERBE/i 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 


1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

167 

25 

0 

0 

1 

0 

0 

0 

0 

0 

3 

0 

0 

557 

602 

0 

0 

5 

72 

0 

0 

0 

3 

4 

0 

0 

22 

354 

0 

0 

1 

241 

0 

0 

2 

161 

4. low 

clouds 



5. 

mid. 

.high 

clouds 

6, 

high clouds 

AVHRR/j 

ERBE/i 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

123 

19 

0 

0 

13 

4 

0 

0 

3 

0 

3 

0 

0 

124 

153 

0 

0 

113 

246 

0 

0 

104 

238 

4 

0 

0 

4 

9 

0 

0 

7 

80 

0 

0 

15 

528 
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Table IV: Averaged standard deviations 




°Mi 

C%] 



Pi 

= 1/Npixel E /100 

Mi 



scene 

pi xel 

-range in 

ffSW 

p sw 

plw 


ERBE 

scan-line (62 pixel ) 

[W/m**2) 

[%] 

m 


12 - 

16 (backward scatter) 

2.1 

2.1 

0.085 

1 






clear 

30 - 

34 (sub - track) 

3.2 

3.9 

0.092 

ocean 

47 - 

51 (forward scatter) 

3.8 

4.3 

0.096 


12 - 

16 

6.3 

4.9 

0.2 

2 

partly 

30 - 

34 

9.0 

8.0 

0.2 

cloudy 

47 - 

51 

9.3 

7.0 

0.2 


12 - 

16 

7.5 

4.1 

0.3 

3 




. 


mostly 

cloudy 

30 - 

34 

15.7 

12.7 

0.4 

47 - 

51 

16.3 

8.1 

0.3 


12 - 

16 

2.1 

0.5 

0.1 

4 






overcast 

30 - 

34 

2.9 

0.6 

0.2 


47 - 

51 

3.0 

0.8 

0.1 




Spectral response of ERBE scanning radiometers channel 
(Smith et al . , 1986) . 



REFERENCE 

LEVEL 

FLUXES 


REGIONAL 

AVERAGES 


T 

M 


2. Inversion processing steps (Smith et al 


. , 1986). 







GEOTYPE: OCEAN 
SCENE ID 


LONGWAVE 

MEASUREMENT 



3. Schematic of maximum likelihood estimate (Smith et al . , 1986) 






Radiative transfer calculations 



4. Schematic flow of simulation approach 











(km) 



RH (%) T (K) 


5. Vertical profiles of temperature and relative humidity for five 
atmospheric models (McClatchey, 1983) used in this study. 
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Spectral responses of NOAA 9 AVHRR channels 1 , 2 and 4 (from 
NOAA Operational Satellite User's Manual). 



Model-filtered SW radiances 
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Computed means arid standard deviations for filtered radiances as 
function of viewing zenith angle for 
a) ERBE shortwave channel and b) longwave channel. 
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N0AA9-ERBE DATA 


Day : 4/4/85 
Colat. : 72.2 - 91.6 
Lon. : 215.1 - 220.1 


G 0 : 32.5 - 37.5 
Q x : 0.0 - 10.0 
4> : 150.0 - 180.0 

• : CLEAR OCEAN 

o : OVERCAST 



Model : TRO 

Scenes 
1 : OCEAN 

5 : LT CLOUD 

6 : MT CLOUD 

7 : HT CLOUD 

13 : LM CLOUD 

14 : MM CLOUD 

15 : HM CLOUD 

19 : MK CLOUD 

20 : HK CLOUD 


8 
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8. Bi spectral distribution of ERBE measurements aboard NOAA 9 spacecraft 
and modeled one-sigma ellipses within for a sample case. 

For April 4, 1985, area 1197 (see text). 




Location of orbit tracks for data periods. 



a 


12.0 ~ 
10.8 - 


NOAA9-AVHRR DATA (GAC) 

Day : 1 4/4/85 
Coiat. : 80.2 - 85.1 



: 226.4 - 227.6 

0 O : 37.5 - 42.5 
0, : 50.0 - 60.0 
4> : 150.0 - 180.0 

Model : TRO 
Scenes 
1 : OCEAN 

5 : LT CLOUD 

6 : MT CLOUD 

7 : HT CLOUD 

13 : LM CLOUD 

14 : MM CLOUD 

15 : HM CLOUD 
• 19 : MK CLOUD 

20 : HK CLOUD 
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N0AA9-AVHRR DATA (GAC) 


Day 



: 14/4/85 
80.3 - 85.3 
• 354.1 - 355.3 

0 O : 37.5 - 42.5 
0, • 50.0 - 60.0 
<P : 150.0 - 180.0 

Model : TRO 

Scenes : 

1 : OCEAN 

2 : LAND 

8 : LT CLOUD 

9 : MT CLOUD 
10 : HT CLOUD 

13 : LM CLOUD 

14 : MM CLOUD 

15 ■ HM CLOUD 

19 : MK CLOUD 

20 : HK CLOUD 
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10. Bispectral distribution of AVHRR measurements aboard NOAA 9 
and modeled one-sigma ellipses within for a sample case. 

a) for April 4, 1985, area 1197 

b) for April 14, 1985, area 1219 (see text). 
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AVHRR VIS — CHANNELS RATIO 
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12. Histogram of ratio Channel 2/Channel 1 for AVHRR data for differing 
scene types. 







14. Histograms of sums of matrix diagonals 
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15. Branching of computational flow due to scene identification. 
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16. Schematic of scene identification probability computation. 
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